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Abstract: The investigation of the scalar gluonium correlator is interesting because it 
carries the quantum numbers of the vacuum and the relevant hadronic current is related 
to the anomalous trace of the QCD energy-momentum tensor in the chiral limit. After 
reviewing the purely perturbative corrections known up to next-next-to-leading order, the 
behaviour of the correlator is studied to all orders by means of the large-/3o approximation. 
Similar to the QCD Adler function, the large-order behaviour is governed by the leading 
ultraviolet renormalon pole. The structure of infrared renormalon poles, being related to 
the operator product expansion are also discussed, as well as a low-energy theorem for the 
correlator that provides a relation to the renormalisation group invariant gluon condensate, 
and the vacuum matrix element of the trace of the QCD energy-momentum tensor. 
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1 Introduction 

The investigation of two-point correlation functions of hadronic currents plays a very im- 
portant role in phenomenological applications of QCD. The field gained much impetus after 
the seminal work of Shifman, Vainshtein and Zakharov [1, 2] in which mesonic correlation 
functions were studied. Gluonic correlation functions were examined in refs. [3, 4], but 
in subsequent years much less work went into the investigation of their phenomenology, 
presumably because the situation on experimental data for the spectra was and still is 
rather difficult. This is in particular so for the scalar gluonium sector, which carries the 
quantum numbers of the vacuum and where strong mixing with scalar meson states is 
expected. A selected set of QCD sum rule analyses includes the works [5-9] where also 
further references can be found. (See also the "Note on scalar mesons" in the Review of 
Particle Physics [10].) 

From the theoretical perspective the investigation of the scalar gluonium correlation 
function has several interesting aspects. First of all, the respective current is given by 
the gluonic piece in the QCD Lagrangian. Next, the corresponding renormalisation group 
invariant (RGI) current constitutes the gluonic contribution to the anomalous trace of the 
QCD energy-momentum tensor, one of the sources of the breaking of conformal symmetry 
[11]. And finally, the expectation value of these currents with respect to the full QCD 
vacuum represents the so-called gluon condensate, one of the central parameters in the 
framework of QCD sum rules. The gluon condensate then also emerges in a low-energy 
theorem for the scalar gluonium correlator [4]. 

The definition of the gluon condensate is plagued with an inherent ambiguity, being 
related to renormalon singularities of the Borel transform of the correlation function in 
the complex Borel plane [12, 13]. The ambiguity arises because a prescription has to be 
chosen of how to treat the infrared (IR) renormalon singularities on the positive real axis. 
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in order to define the Borel integral. On the other hand the IR renormalon singularities 
are connected to the power corrections in the operator product expansion (OPE) and 
their general structure can be analysed with the help of the renormalisation group [13-18]. 
Still, complete results for the Borel transform of two-point correlators are most often only 
available in simplified models, like the large-nj or relatedly the large- /3o approximation, 
where Uf is the number of fermion flavours and /3o the leading coefficient of the QCD f3- 
function. Nonetheless, valuable information on the structure of the correlators and their 
large-order behaviour can be gleaned from such approximations. In the following work, 
therefore, the above aspects shall be examined for the scalar gluonium correlation function. 

To begin, in section 2 the general structure of the perturbative series for the scalar 
gluonium correlator will be discussed. To this end, two choices for the hadronic current 
will be compared. The first consists of the gluonic term in the QCD Lagrangian. However, 
as this current is not renormalisation group invariant, the correlator cannot directly be 
connected to an observable quantity. For this reason then a second, related correlator 
will be investigated, which is constructed from the RGI combination that appears in the 
trace of the QCD energy-momentum tensor. The spectral function (imaginary part) of this 
correlator is an observable and also another physical quantity can be constructed which 
bears similarities to the Adler function for the mesonic vector correlator. The analytically 
available information on the perturbative coefficients of these correlators will be reviewed. 

In section 3, the large- /3o approximation for the scalar gluonium correlator will be com- 
puted and discussed. In particular the structure of the ultraviolet and infrared renormalon 
poles will be investigated, which are related to the large-order behaviour of the perturba- 
tive expansion of the correlation function and the structure of higher-dimensional operator 
corrections in the OPE. Furthermore, based on RG arguments, the structure of the leading 
IR renormalon singularity will be derived beyond the large-/3o approximation. Section 4 
finally discusses a low-energy theorem for the scalar gluonium correlator [4] and its relation 
to the gluon condensate, before the work is concluded in section 5. 

2 Scalar gluonium correlator 

The most basic two-point correlation function that is relevant for the study of scalar gluo- 
nium can be defined as 

^G^iQ^) ^ i J dxe*-?^ {n\T{JG{x)JG{0)m , (2-1) 

where the gluonic current is given by Jg{x) = G^^{x) G^"°'{x), G^^{x) is the QCD field- 
strength tensor and \Q) represents the full non-perturbative QCD vacuum.^ 

Analogous to the mesonic vector correlator (see e.g. ref. [18]), the purely perturbative 
part of 11(^2(5^) has the following general structure: 

n^I(.) = - ± a; g L^ ^ - In ^ , (2.2) 

n=0 k=0 ^ 

^Operator products at the same space-time point are assumed to be renormalised by (modified) minimal 
subtraction in dimensional regularisation. 
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with = a{fj?) = as{fj,)/'ir and fi the renormahsation scale. 11(^2 (s) itself is not a physical 
quantity. However, after multiplying an appropriate function in the coupling that is related 
to the QCD /3-function and will be discussed in detail below, the spectral function is 
physical, as well as the third derivative of 11^2(5) with respect to s, because for large s, 
n^J(s) scales as and three subtractions are required. Let us first concentrate on the 
spectral function pq2{s), related to IIq2{s), 

Pq2{s) = -lmnG2{s + iO) . (2.3) 
vr 

Up to second order in the explicit expansion of the perturbative spectral function 
p^{s) is found to be 

Pc^is) = 47^2 ^ ^[go.i + (ci.i + 2ci.2-^)«M + Ki-^^C2,3 + 2c2,2^ + 3c2,3-^^)a^J , (2.4) 

where s > and L = ln(s//i^). Like for the Adler function, the perturbative coefficients 
Cnfi are unphysical constants that have to drop out for every physical quantity, the c„^i 
can be chosen as the only independent coefficients, and, as will become evident from what 
follows, all other Cn,k for k > 1 can be determined from the renormahsation group equation 
(RGE). Comparing eq. (2.4) with the analytic computation of eq. (10) of ref. [19] (for an 
earlier next-to-leading order calculation see ref. [20]), and setting Nc = 3, the first three 
independent coefficients are found to be: 

^—1 ^— 73 7^ ^ — 37631 495 /- /7189 5 a \^ , 127 ^2 /r, rN 

In order to deduce constraints from the RGE, we have to work with a renormahsation 
group invariant current. In the chiral limit, where the operator Jg{x) does not mix with 
ruq q{x) q{x) or m^, such a current can be chosen to be [21] 



Jg(x) = ^Jg{x) = ^G;,{x)G^^%x) « aG;,{x)G^^^{x). (2.6) 



Here, the QCD /3-function is defined as 

/3(a) ^ -p^ = /3ia2 + /32a3 + ... . (2.7) 
ap 

Explicitely, at Nc = 3, the first coefficient of the /^-function reads /3i = 11/2 — nj/3. Let 
us remark that in the chiral limit Jg{x) is proportional to the trace of the QCD energy- 
momentum tensor 9li{x) [21-23], with the relation being given by 

In analogy to nfj2 (g^) of eq. (2.1), we can define the two-point correlator for the current 
Jg{x), which expressed in terms of (g^) takes the form: 



^G<i') = (^) ^G^ir). (2.9) 
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The spectral function Pc^is) = lmIlQ2{s + i0)/7r corresponding to this correlator now is a 
RGI "physical" quantity, which satisfies a homogeneous RGE. Working out the constraints 
from the RGE, up to next-next-to-leading the following relations for the coefficients c„^fc 
are found: 



Cl,2 



C2,2 



/32 Co,l - ^^ Ci,i , 



/3f 

C2,3 = X^O,! 



(2.10) 



From the logarithmic terms of eq. (10) of ref. [19] it can easily be verified that these relations 
are satisfied. 

The results given above can also be compared to corresponding expressions of ref. [24]. 
Taking into account the fact that the definition of the gluonic current in [24] is Jq^^{x) = 
nJcix), and setting the number of quark flavours nj = 3, one obtains 

^ _ o „2 r 1 I 659 ^ I / 822569 465 a \„21 

ao--2a [1 + "36- a + (^592 ^Qs) a \ , 



ai= 2a=^[| + 



2105 
32 



0-2 



(2.11) 



with ao, oi and 02 having been defined in eq. (30) of ref. [24]. The above expressions are 
found in agreement to eq. (31) of [24], there given partially numerically. 

As has been noted before, since IIq2 (s) behaves proportional to for large s, it satisfies 
a dispersion relation with three subtraction constants. 



Ug2 (s) = Ug2 (0) + s (0) + - (0) + 



Paris') 



iO) 



ds'. 



(2.12) 



where the prime denotes a derivative with respect to s. Therefore, besides the spectral 
function pQ2{s), taking three derivatives of IIq2(s) another function Dq2{s) which is a 
physical quantity can be defined in analogy to the QCD Adler function: 



Dg2{s) 



d'nc,2(s) 

ds3 



(2.13) 



Employing eq. (2.2), the general structure of the perturbative expansion for Dq2 (s) is found 
to be 



1) 



_ ^ 

+ 3{k - 



2 00 



n+l 



2L 



n=0 k=l 



l){k - 2)L 



k~3 



(2.14) 



Resumming the logarithms in D^{s) with the scale choice pp' = — s, the general expansion 
as well as the perturbative expression up to second order are found to be 

Dins) = (1^)' E [^-^ + 3^-^ + 3^-3] ai-sT 

4 



-^a2(-s) 



vr 



n=0 



161 , , 
l + ^«(-s) + 



/ 1269361 495 « 2/ ^ 



5808 



, (2.15) 
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where in the second hne the exphcit expression for nj = is provided. One notices that 
the perturbative higher-order Us corrections are large. To obtain a reasonable behaviour 
of this perturbative series, one should have as{—s) as small as roughly 0.12 which requires 
scales as large as about the Z-boson mass scale Mz, and thus precludes a low-energy sum 
rule analysis. Nonetheless, one may proceed to investigate the general properties of the 
gluonium correlation function. 

Since aspects of the perturbative series related to the renormalon ambiguity of the 
gluon condensate shall be studied, it is also of interest to consider the large-/3o approx- 
imation (or large-/?! approximation in my notation, though I shall keep the historical 
terminology). In this approximation the Adler-like function of eq. (2.13) turns out to be: 



Numerically, the large- /3o approximation is reasonably close to the full result, 11 versus 
14.6 at order and 104.6 versus 144.2 at order a^, so that it appears to make sense to 
work out the complete large- /3o approximation to D^{s). The behaviour of the series in 
large-/3o is slightly better than in full QCD as the perturbative coefficients are somewhat 
smaller. The large-/3o approximation to all perturbative orders shall be developed in the 
next section. 

3 Large-/3o approximation 

Like in the case of mesonic correlation functions, the large-/3o approximation can be cal- 
culated by first computing fermion-chain diagrams that correspond to the limit of a large 
number of fermion flavours nj and then applying the so-called naive non-abelianisation, 
that is replacing the nj dependence by the leading coefficient of the QCD /3- function, which 
then also incorporates a gauge invariant set of gluon-loop diagrams [25, 26]. In practice, 
one may either explicitly calculate the required fermion-chain diagrams for 11^2 or proceed 
directly to a computation of the Borel transform of IIq2 or Dq2 . As below the structure 
of the Borel transform B[Dq2{s)]{u) will be discussed in relation to the OPE, first of all 
this Borel transform shall be presented. A discussion of the large-/3o result for 11^2(5) is 
relegated to appendix A. 

Technically, the computation of the large- /3o approximation for the gluonium correlator 
proceeds through inserting chains of fermion loops in both gluon propagators of the lowest 
order gluon loop diagram. Employing the Dyson-resummed form of the gluon propagator, 
together with the convolution theorem for the Borel transform (see e.g. eq. (2.8) of ref. [27]), 
one arrives at 



ia2(-s)[l + 2/3ia(-s) + i/32a2(- >.) + ...] . (2.16) 



B[Dg2{s)]{u) 



( 




—II 



r(i + n) 
r(4 - u) 



u 



I 



du[2 — u + u{u — u 



r{2 -u)r {2 -u + u) 

r{l + u)T{l + u-u) 
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fli 


02 


as 


a4 


05 


ae 


1 


2 


2.30556 


1.76444 


1.14045 


0.55870 


0.30433 


07 


as 


ag 


aio 


an 


ai2 


ais 


0.09359 


0.07845 


-0.00719 


0.03581 


- 0.02458 


0.02891 


-0.02726 


ai4 


ai5 


ai6 


ai7 


ai8 


ai9 


(120 


0.02789 


- 0.02764 


0.02774 


- 0.02770 


0.02772 


- 0.02771 


0.02771 



Table 1. Numerical expansion coefBcients a„ of eq. (3.3) up to the 21st order. 



where C is a scheme-dependent constant which in the MS scheme takes the value —5/3. 
Assuming a positive coupling ag, the relation between the Borel transform B[Dq2]{u) and 
the Adler-like function Dq2 is given by 



00 



Dc^a) = —jdue B[Dg2]{u) . (3.2) 





Without further specification, the Borel-integral on the right-hand side is only well defined 
if B[Dq2]{u) has no poles or cuts on the positive real axis. Similarly to the case of the QCD 
Adler function, as will be discussed in more detail further below, this is not the case for the 
expression of eq. (3.1). Still, the perturbative expansion, corresponding to an expansion of 
B[Dq2]{u) in powers of u, is well defined, and takes the form 

S[Z?G.(.)](..) = ^-£3-^ ^a„n"+i 



n=0 



7r3/3i 



u 



+ 2n2 + in^ + (lf-|)n^ + ...l , (3.3) 



where in the second line the first few orders have been given explicitly. Higher-order 
expansion coefficients a^ are collected in table 1 in numerical form up to the 21st order. 
Together with the formula 

duu^'e'^ = n! ^ , (3.4) 







2 



leading to the perturbative expansion 

^ia.ge-,o(,) = ^^a\-s) f; ^^a.[/3la(-.)]^ (3.5) 

n=0 

it is a straightforward matter to recover eq. (2.16) for DQ2{a) in the large-/3o limit. Eq. (3.5) 
also clearly displays the factorially divergent asymptotic behaviour of the perturbative 
series. 
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Next, the singularity structure in the Borel plane, also known as renormalons, shall 
be investigated. Via the factor r(l + u), B[Dq2]{u) has ultraviolet (UV) renormalons for 
all negative integer u. The UV renormalon singularity at u = — 1 is also the closest to 
the origin and thus dominates the perturbative series at large orders. In the vicinity of 
u = —1, the behaviour of the Borel transform reads 

(jv^ -^\ J93_ c 
B[Da.]{n) "^-^ - ff^ , (3.6) 

vr-'pi (1 + u) 

from which the perturbative coefficients at large orders can be deduced. In the MS scheme, 
dominance of the leading UV renormalon singularity sets in at about the 12th order, and 
inspecting table 1, it is found to be in accord with the behaviour of the coefficients at 
sufficiently large n. 

The singularities at positive n, the infrared (IR) renormalons, are connected to the 
power corrections in the OPE. For physical correlation functions they are generally expected 
to start at ti = 2 which corresponds to the dimension-4 gluon condensate. In the case at 
hand, however, the correlator 11(^2 (s) itself is of dimension four, and therefore the gluon 
condensate only appears as a subtraction constant which vanishes for the physical correlator 
Dq2. This is also seen from the fact that the Borel transform B[Dq2]{u) stays finite at 
u = 2, and will be further discussed in connection to a low-energy theorem for IIq2 in the 
next section. Hence, potential IR renormalon singularities only start at integer u > 3. 

Even though from RG arguments generally power-like IR renormalon singularities are 
expected [13-18], in large-/5o for B[Dq2]{u) they turn out to be only logarithmic. To 
understand the origin of this behaviour, first the generic case shall be investigated. The 
derivation is completely analogous to the one for example presented in section 5 of ref. [18], 
the only difference being that now the correlator has an additional multiplicative factor 
of af , with 6 = 2 for Dq2. The structure of an IR renormalon pole corresponding to an 
operator of dimension d is then found to be 

B[dM^) = j^-^[l + Oip-u)], (3.7) 

where are normalisations of the renormalon poles that cannot be determined from the 
renormalisation group, and 

Here, is the leading order anomalous dimension of the operator O^. 

At It = 3, the only contributing operator corresponds to the dimension-6 gluonic 
condensate (0|g'/"^'^G^^G'^'^G'^^''^|il), whose leading anomalous dimension was calculated 
to be IgQ-i = {7Nc + 2nj)/6 in [21, 28]. Hence, in the absence of quarks, that is nj = 0, 
the exponent in the denominator of eq. (3.7) turns out to be 1 + 7 = 108/121 (independent 
of Nc). In the large-/3o limit, on the other hand, the /32-term vanishes and "y^^ls/Pi = — 1. 
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Thus, the exponent 1 + 7 is found identicahy zero. This is reflected in the fact that close 
to n = 3 the Borel transform (3.1) behaves as 

B[DaA{u) - ^^5^ 3e-3^ In [l - |] , (3.9) 

only displaying a logarithmic singularity. Such a logarithmic behaviour of the leading IR 
renormalon was also observed in ref. [17]. A further consequence of the softening of the 
leading IR renormalon singularity is the fact that at no order it dominates the perturbative 
coefficients. This is in contrast to the QCD Adler function for which in large-/3o the leading 
IR singularity at n = 2 gives the most important contribution at intermediate orders before 
the large-order behaviour governed by the UV singularity at u = —1 sets in [13].^ 

At positive integer p > 4, the Borel transform B[Dq2\{u) again behaves logarithmically 
like ln[l — u/p], but this time, due to the factor l/r(4 — u) in eq. (3.1), these logarithms 
are further suppressed by an additional factor {p — u). Thus, at positive integer n > 4, 
the Borel transform does not posses further singularities, but just the appearance of new 
cuts which reflect the presence of higher-dimensional contributions in the OPE. In order to 
properly define the Borel integral of eq. (3.2), hence a prescription should be chosen which 
circumvents the singularity at n = 3 and the cuts at all integer n > 3. This prescription 
then also gives a meaning to the QCD condensate parameters. 

4 Low-energy theorem 

Even though IIq2{s), and therefore also 11^2(0), are not physical quantities, the question 
arises if a closed expression can be given for the subtraction constant 11(^2(0) of the dis- 
persion relation (2.12), similar to the case of the pseudoscalar correlation function ^5(5) 
to be discussed further below. Such relations are termed low-energy theorems (LET) and 
they are interesting as they provide additional constraints on the theory in question. In 
fact, already in the seminal work [4], a LET for 11(^2 (s) was derived which reads: 

= lim %.{s) = ^ Warn . (4.1) 

The LET can also be expressed in terms of the correlator I\.e{s) for the trace of the energy- 
momentum tensor Oli{x) of eq. (2.8), in which case it assumes the particularly simple form 

n,(0) = hm Heis) = -4(0|e^|0) . (4.2) 

The gluon condensate which appears on the rhs of (4.1) and (4.2) has a renormalon am- 
biguity. This is for example seen in the large- /3o approximation for the Borel transform 
of 11(52(5) which, in contrast to B[Dq2\{u), has a UV renormalon pole at n = 2. (See 
also appendix A.) This pole provides an ambiguity of order Aqqj-, which cancels against a 
corresponding ambiguity in the definition of the gluon condensate. Nevertheless, this does 
not imply that the gluon condensate is an observable, since the subtraction constant is not 
directly related to the spectrum of ng.2(s). 

^The latter behaviour is also expected for the Adler function in full QCD [18]. 
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The situation is analogous to the case of the pseudoscalar correlator ^'5(5). Here the 
subtraction constant \I'5(0) is related to the quark condensate {^l\qq\Q), and the low-energy 
theorem can also be derived for a hadronic representation of the currents which yields the 
Gell-Mann-Oakes-Renner (GMOR) relation [29]. The derivation of the GMOR relation can 
be generalised in the framework of chiral perturbation theory (xPT) [30], and the fact that 
the subtraction constant ^'5(0) is not a physical observable then manifests itself through 
the appearance of the unphysical coupling H2 at higher orders in the xPT expansion 
[31], reflecting the dependence on the short-distance renormalisation. 

It might be perhaps interesting to study a GMOR-like relation also for the scalar 
gluonium correlator IIq2{s). The framework for such an analysis is available from the work 
of Donoghue and Leutwyler [23] , who derived chiral expansions for correlation functions of 
the energy- momentum tensor. As the spectrum of the scalar sector is not very well known, 
and a strong mixing between gluonic and mesonic states is present, the analysis is however 
expected to be more involved and will be left for the future. Still, as for example the 
recent dispersive analysis [32] of the scalar / = sector demonstrates, the situation is not 
completely hopeless. Additionally, the LET may open the route to a lattice determination 
of the gluon condensate. 

5 Conclusions 

The investigation of the scalar gluonium correlation function has several interesting theo- 
retical aspects. The relevant interpolating current represents the gluon term in the QCD 
Lagrangian and its renormalisation group invariant counterpart is proportional to the glu- 
onic piece in the QCD energy-momentum tensor. Furthermore, a low-energy theorem for 
the correlator provides a relation to the gluon condensate, an important parameter in the 
framework of QCD sum rules. 

The general structure of the perturbative expansion for the scalar gluonium correlator 
was discussed in section 2, and explicit analytical results are available up to next-next-to- 
leading order. Furthermore, the physical correlation function Dq2{s) was defined, which 
is the analog of the QCD Adler function in the case of the vector current correlator. It 
turns out that the perturbative corrections for Dq2 (s) are rather large, for which reason a 
phenomenological analysis of this correlator would appear questionable, leaving aside the 
additional complication of the sparse knowledge on the physical spectrum. 

Still, from the theoretical perspective a general analysis of the behaviour of the per- 
turbative expansion, also at large orders is interesting, because on the one hand it sheds 
light on the asymptotic nature of the perturbative series, and on the other hand the sin- 
gularity structure of the Borel transform of Dq2{s) in the positive Borel plane is related 
to higher-dimensional operator corrections in the operator product expansion. 

The all-order perturbative behaviour of current correlators can be studied in the large- 
nj, or relatedly, the large- /3o approximation. To this end in section 3 the Borel transform 
B[D(j2]{u) was calculated, from which it is a simple matter to derive the large-/3o approxi- 
mation of the correlator Dq2 (s). Like in the case of the QCD Adler function, the large-order 
behaviour is governed by the leading ultraviolet renormalon singularity at u = —1, and in 
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the MS scheme dominance of this contribution takes over at about the 12th perturbative 
order. 

Contrary to the Adler function, Dq2{s) does not possess an infrared renormalon pole 
at u = 2, because the gluon condensate only appears in the subtraction constant 11(72(0), 
which due to the derivatives vanishes for the physical correlator Dq2{s). This statement 
should also hold true in full QCD. Furthermore, in the large-/3o approximation the leading 
infrared renormalon singularity at u = 3 is only logarithmic, and it is observed that at no 
order this renormalon pole provides the dominant contribution to a certain perturbative 
coefficient. The divergent structure is, however, different in full QCD, for which the power- 
like renormalon singularity at u = 3 was explicitly derived from renormalisation group 
considerations. 

An interesting limiting case are the low energy theorems of eqs. (4.1) and (4.2), since 
they are solely expressed in terms of the renormalisation group invariant gluon condensate. 
This does, however, not imply that the vacuum expectation value of the (even RGI) gluonic 
operator is a physical quantity. Like in the case of the pseudoscalar correlation function, 
the subtraction constant is not directly related to the physical spectrum and is expected 
to depend on the short-distance renormalisation. Still, the low-energy theorems may open 
additional routes to extract a phenomenological value for the gluon condensate. 

These considerations bear some relevance for the question of the QCD contribution 
to the vacuum energy, as such a contribution is given by the vacuum expectation value of 
the trace of the energy-momentum tensor. From the discussion above it appears as if at 
least in pure QCD this contribution is unphysical, as it depends on the renormalisation 
prescription. An immediate interesting question arises: what happens if QCD is coupled to 
gravity? As the expectation value of the energy-momentum tensor between particle states 
is physical, and a gravitational field can change particle number, then also its vacuum 
expectation value should become physical. The investigation of this question is beyond the 
scope of the present article, but should be pursued in the future. 
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Appendix A: Large-/3o approximation for 11(52 (s) 

In the following, the large-/3o approximation for the correlator IIq2{s) of eq. (2.2) will be 
provided and discussed. The expression which arises from the explicit computation of the 
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double bubble-chain diagrams in a general space-time dimension D = 4 — 2e reads 



ngr"*(s) 



m n 

■Y.T. 

k=0 /=0 



EE 

^ m=0 n=0 



2e 



m+n 



m\ n 



■6 5(2 - e, 2 - e)r(l -Fe) 



47r^ 



2\e 



2 - {?, + k + l)£ + {k + + l)e'^ Gki{e) 



(A.l) 



where l/e = 1/ £ — + ln(47r), B{x, y) is Euler's beta function, and the function Gkiie) is 
found to be: 



Gki{e) 



r[2-(fc + i)g]r[2-(/ + i)g] 

r[l + A;e]r[l + /e]r[4- (A; + / + 2)e] 



r[-2 + {k + I + l)e] 



(A.2) 



The independent coefficients c^'j in the large-/3o approximation can then be obtained by 
working out the terms linear in the logarithm L = ln{—s/fj?) from eq. (A.l). 

In order to be able to compare to the large-/3o approximation for the function Dq2 (s) 
of eq. (3.5), one also requires the coefficients cj^2 c^fg, that is the terms quadratic and 
cubic in L. These can either be calculated from a direct expansion of (A.l), or through 
the following relations, valid in the large-/3o approximation, which are derivable from the 
RGE: 



.'/Jo 



1 

12 



4 

n + l 



n-1,1 ' 



+ 1/3^^2,1, 



(A.3) 



(A.4) 



where [r] denotes the integer part of r. Employing all results, one can verify that the 
relation 

{n + iy 



2" 



Wo 

n,3 



(A.5) 



is satisfied, which provides a good check of the computation. 

A final remark concerns the Borel transform of (s). Through a factor T[u — 2 + £), 
at -D = 4 it contains UV renormalon singularities for all integer u < 2. The ambiguity 
inflicted by the rightmost singularity at u = 2 corresponds to the corresponding one of 
the gluon condensate in the subtraction constant 11(^2(0). The UV renormalon for the 
gluon condensate at n = 2 has also been obtained in ref. [33]. Considering in addition the 
momentum dependence, and using the relation 



- s-^s2-"-^r(u-2 + e) 



ds3 



'T{l+U + £) 



(A.6) 



it becomes clear that at D = 4 the Borel transform of the physical correlator Dq2 (s) only 
has UV renormalon singularities for u < — 1, as discussed in section 3. 
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